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' I. INTRODUCTION. 

CNi : 

^ One of the most important problems of theoretical physics in the last years was to answer the question: How and 

1 ^ in what circumstances a quantum system becomes classical ? [0 . In spite of the great effort made by the physicists 

to find the answer, the problem is still alive |^ and we are far from a complete understanding of many of its most 
fundamental features. In fact the most developed and sophisticated theory on the subject, histories decoherence, is 
I not free of strong criticisms Q . 
^ . For conceptual reason we will decompose the limit quantum mechanics classical mechanics in two processes: 
^SJ ' quantum mechanics — > classical statistical mechanics and quantum statistical mechanics — > classical mechanics. There 
T-H , is an almost unanimous opinion that the first process is produced by two phenomena: 

' i. Decoherence, that in quantum systems, restores the boolean statistic typical of quantum mechanics and 
I ii.- The limit h ^ 0, that circumvents the uncertainty relation at the macroscopic level and allows to find the classical 
3^ ' behaved density functions via de Wigner integral. 

I The second process is produced by these phenomena plus the localization (or production of correlation) phenomenon 
— ... , 1^. We will discuss at large the first process and briefly deal with the second at the end of each section. 
^ 1' The techniques to deal with the two phenomena related with the first process are not yet completely developed. 
I , One of the main problems is to find a proper and unambiguous definition of the, so called, pointer basis, where, 
decoherence takes place. 

rr Our contribution to solve this problem is based in old ideas of Segal (_|^, Q) and van Howe reformulated by 

K>" Antoniou et al. We have developed these ideas in papers |^ and jl^ where we have shown how the Riemann- 
KJi Lebesgue theorem can be used to prove the destructive interference of the off-diagonal terms of the state density 
matrix yielding decoherence. Using this technique we have found decoherence and the classical statistical equilibrium 
' limit in simple quantum systems Q| where we have defined the final pointer basis in an unambiguous way. g. Also 
the localization phenomenon appears in some cases. 

This paper is devoted to give two examples of the method introduced in papers ||^ and , that we will briefly 
review in section II, and used it to find a general solution for the quantum-classical limit in paper 

In section HI, the method will be used to solve the problem known as "the Mott problem" after the name of Sir 
Neville Mott, probably the first one who studied the subject Let us consider a radioactive nucleus, placed at the 
origin of coordinates O, inside a bubble chamber. We will see that classical radial trajectories appear in the chamber 
due to the emitted out-going particles. We must explain this phenomenon. Theoretically we have a timeless structure, 
since the wave function, iplx) satisfies the eigenvalue equation: 

i/V = wV' (1) 

where: 



'^The relation of our method with the histories decoherence is studied in paper Q. They turn out to be equivalent, but in our 
method the final pointer basis is more properly defined. 
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H = - — -—r^— + ^— + W(r) (2) 

being W{r) the spherically symmetric potential barrier, which is the external wall of a potential well, and such that 
limr^oo W{r) = 0. This will be our model for the nuclear forces Let us observe that there is no trace of the time 
in these equations. 

The main problem is that, even if the symmetry is a spherical one, there is no a priori reason to explain why the 
classical trajectories are radial. Really we have the following list of facts to explain: 

i. -Why the nucleus quantum regime becomes the classical regime for the classical trajectories 

ii. - Why these classical trajectories are radial. 

iii. - Why the notion of time, necessary to explain the motion of the classical radial particles, appears in a timeless 
formalism, i. e. in the one of eq. (|l]). 

iv. -Why there are only outgoing motions 0. 

At this stage we must observe that this set of problems is very similar to the one of Quantum Cosmology, where, 
in fact, we must explain the outcome of the classical regime, the appearance of time, the nature of the classical 
trajectories in superspace, and the direction of the corresponding motion: i. e. the arrow of time. Then several of 
the most important quantum universe problems, that we will discuss in the second example, are already contained in 
our humble Mott model. 

So, in section IV, we will consider our second example: the quantum cosmology problem, since the appearance of 
a classical universe in quantum gravity models is the cosmological version of the first problem. Then, decoherence 
must also appear in the universe JlSt . 

In this paper, using our method, we will solve the two examples and we will find: 

i. -Decoherence in all the dynamical variables and in a well defined final pointer basis. 

ii. -A final classical equilibrium limit, when ?i — > 0, in such a way that the Wigner function F^' = pi"^^ of the 
asymptotic diagonal matrix can be expanded as (e. g., in the cosmological problem): 

pl^')([x], [k]) = J P{0[a]p53a„](W> H)rf{/}d[a] (3) 

where pj-'^j^j^^j ( [x] , [k]) is a classical density strongly peaked ^ in a trajectory defined by the initial coordinates a and 
the momenta I and P{i}[a] is the probability of each trajectory. As essentially the limit of quantum mechanics is not 
classical mechanics but classical statistical mechanics this is our final result: the density matrix is translated in a 
classical density, via a Wigner function, and it is decomposed as a sum of densities peaked around all possible classical 
trajectories, each one of these densities weighted by their own probability ^ 

Thus our quantum density matrix behaves in its classical limit as a statistical distribution among a set of classical 
trajectories. Similar results, for the cosmological case, are obtained in papers and jist . 

We will reelaborate these conclusions in section V. 



II. REVIEW OF THE METHOD. 



In order to go from the quantum to the classical statistical regime two new properties must appear: 
i. -Decoherence: The density matrices, that contain quantum interference terms, must become diagonal, in such 
a way that these interferences would be suppressed. Then the quantum way to find probabilities of exhaustive events 
(i. e.: adding the corresponding amplitude and computing the norm) become the classical way: just adding the 
probabilities. 



^But let us complete the first example discussing the role of the global atmosphere of the bubble chamber in our problem. 
This role is non essential and really just incidental. In fact, if there were no bubble chamber, at least for r —> oo the classical 
trajectories would also be radial, since in this case, any small detector located far enough from the origin would find radial 
motions. Furthermore it would be very difficult to believe that these radial motions are produced only by the small detector. 
Then the radial motion exists even if there is no bubble chamber and it is a mistake to consider that it is the bubble chamber 
that, acting as an environment, produces the radial structure. But we must remark that there is a measurement in both cases 
and therefore the measurement processes is essential. 

^Precisely: peaked as allowed by the uncertainty principle. 

^After classical statistical regime is reached correlations will eventually produce a pure classical regime. 
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ii.-The limit h ^ 0: The positions and the momenta (or more generally canonically conjugated dynamical vari- 
ables) can be defined as allowed by the uncertainty principle, but big scales (i. e. when h 0) allow us to consider 
both the position and the momentum as independent dynamical variables, like in classical mechanics. Of course this 
independence is the essential property to find a classical behavior. 

These two closely related properties introduce the classical statistical behavior in the quantum formalism. Let us 
begin with decoherence. 

We are only interested in the scattering states with continuous spectrum, e. g., for the first example, the Mott 
problem, the radial outgoing particles are described by these states. To obtain the hamiltonian eigenbasis of the 
Hilbert space Ti. we can consider e. g. the hamiltonian (^) and construct Lippmann-Schwinger basis [ p7[ , 
(really |a;,/,n-|-)) where to is partially discrete and partially continuous, e. g. it will have some u/q for the ground 
state and a continuous to for the scattering states, but for the moment we will consider only the continuous index uj, 
(since we are only interested in these states |^. Of course we can as well use {\lu— >}■ Then the hamiltonian can be 
diagonalized as: 



H= uj\lu+ X LU + \duj (4) 
Jo 

From this expression we can deduce that the most general observable that we can consider reads: 

/•oc /'OC pOC 

= / 0^\uj+ >< oj + \duj + / 0^^'\uj+ >< uj' + \dujduj' (5) 







where the functions 0^,0^^,^' are regular, namely, the most general observable must have a singular component (the 
first term of the r.h.s. of the last equation) and a regular part (the second term). If the singular term would be 
missing the hamiltonian would not belong to the space of the chosen observables |^ . Then O G O d Ti. (B {Ti. iSi Ti.) 
since some extra conditions must be added This space has the basis {|i^), |tj,ti;')}: 

|w) = ><uj + \, \u;,u;') = >< lu' + \ (6) 

The regular quantum state p are measured by the observables just defined computing the mean values of these 
observable in the quantum states < O >p= Tr{pO) p6| . These mean values can be considered as linear functionals 
p on the vectors O. Then the notion of state can be generalized to any linear functional over O that we can call 
{p\0) In this way not only regular states but singular states can be defined. Moreover, asp must be normalized, 
selfadjoint and positive definite, p G 5 C O' , where 5 is a convex set contained in O' 0, [Q. The basis of O' is 
{(a;|, (cj,cj'|}. These states are defined as functionals by the equations: 

(lu\lu') ~ S{oj — lu'), {lu,lu"\lu',lu"') ~ S{lu — Lu')d{Lu" ~ oj'") (7) 

Therefore a generic quantum state reads: 



oc poo 



P = / Puj{Lu\duj + / / pi^^>{Lu,Lu'\dujduj' (8) 
Jo Jo Jo 

where p^ > 0, p^^c^' = pZ'uj- The states such that p^ ^ p^^^ will be called generalized states Q and those such that 
Puj — Pljlj are the usual regular mixed or eventually pure states. To continue, even if the time is not strictly defined 
(since we have only the eq. (|l|)) let us postulate that there is a symmetry in the system with a symmetry group e~'^*. 
At this level of our reasoning this is a global fact imposed by the structure of the universe where we suppose the 
model is immersed (we will come back to this problem in subsection IIICl, of course this will also be the case in the 
cosmological example, because the problem of time definition is the same). Then the time evolution of the quantum 
state p reads: 



oc roc 



p{t)^ p^{Lu\dLU+ / p^^/e'^'^-'^'^'iLU^Lu'ldLudLu' (9) 

Jo Jo Jo 

As in the statistical level we are considering we only measure mean values of observables in quantum states, i. e.: 



^Bound states are considered in papers Q and [|lC 
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< O {p{t)\0) = 



OO POO 



p^0^duj+ / / p^^,0^^,e-'('^-'^'^'dujdLu' (10) 

"'0 "'0 

using the Riemann-Lebesgue theorem we obtam the weak hmit, for all O G O, p G S: 

lim < O >p(t)==< O >p, (11) 
where we have introduced the diagonal equilibrium state: 

/>oo 

p* ^ Pujii^lduj (12) 
Jo 

Therefore, in a weak sense we have: 

W lim p{t) = (13) 

t — ^oo 

Thus, any quantum state goes to a equilibrium diagonal state weakly, and that will be the result, if we observe and 
measure the system evolution with any possible observable of space O , e. g.: for the first example, with the global 
bubble chamber (or with any one of the local small detectors of the footnote of the introduction). Then, from the 
observational point of view we have decoherence of the energy levels, even if, from the strong limit point of view the 
off-diagonal terms never vanish, they just oscillate (eq. (|9|)). So, from now on, as we will consider the matrix p{t) 
for big e.g.: for the first example, far away from the nucleus, the relevant state will be p*, a diagonal state in the 
energy. 

Some observations are in order: 

i. - The real existence of the two singular parts introduced above is assured by the nature of the problem. The 
singular part of the observables is just a necessary generalization of the singular part of the hamiltonian, which is 
completely singular (eq. (||)). The singular part of the states is the final state for decoherence (eq. (p^). {p\0) is 
just the natural generalization to the continuous of the trace of the product of two finite dimensional matrices (eq. 
©)• 

ii. - To fully understand the phenomenon it is necessary to use generalized states. Therefore, any explanation only 
based on pure wave function states or mixed states is incomplete q 

Having established the decoherence in the energy we must consider the decoherence in the other dynamical variables. 
Before going to the model let us study the general case (we will repeat this reasoning in sections IIB and IVB, with the 
notation corresponding to each case). The diagonal singular component of the eq. (^ (which is equal to p*) is time 
independent, therefore it is impossible that a different decoherence process takes place in this component to eliminate 
the off-diagonal terms, in the other dynamical variables. Therefore, the only thing to do is to try to find if there is a 
basis where these diagonal terms vanish at any time and therefore there is a perfect and complete decoherence. For 
t — > OO this basis in fact exists and it is known as the final pointer basis. 

Let {H,Oi, ...On} be the usual complete set of commuting observables (CSCO) that we are using to make our 
calculations and {\tu,mi, ...,mN+ >} (that we will simply call {jw, toi, mAr >} from now on) the corresponding 
eigenbasis. Then introducing the new indices in eq. ( [l^ ) the equilibrium diagonal state reads: 

PmL.,m„,mi,...,r<('^'™l'-''^A''"^l'--'"4Mw (14) 

From what we have said under eq. (|§) it is: 



Now an incidental question for the first example would be: which energies? In the first example the potential well, surrounded 
by the barrier, W{r) originates unstable levels with energy uj„ which decay with a decaying time 7^^. Inside the well there are 
oscillating waves that can be used to fulfill the boundary conditions at r = 0. When these waves arrive to the barrier they are 
par tially reflected and transmitted The states that tunnel the barrier appear with an energy that peaks strongly at uj„ Q, 
[ |17[ . Therefore the energy of the outgoing particles are energy packets ^ 5{lu — cun), and therefore could be essentially labelled 
with n. But taking into account that the energy spectrum is really a continuous one we will always refer to it as lu. 
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(p^"^ , , )* = P^"^ ' (15) 



therefore this matrix can be diagonahzed and there is a basis {(oj, /i, /a^|} where the matrix reads: 



Now we can define the observables: 



E Ptli>^^^^-Md^ (17) 

and the CSCO {i?, Pi, Pat}, where the singular component is diagonal in the dynamical variables corresponding 
to the observables Pi, Pm from the very beginning. This is the final pointer basis where there is perfect decoherence. 
Final pointer basis is therefore defined by the dynamics of the model and by the quantum state considered, in complete 
agreement with the literature on the subject. The classical statistical limit will be complete if we transform all these 
equations via a Wigner integral as we will in each example using the corresponding notation. 



III. THE MOTT PROBLEM. 
A. Decoherence. 

After these general considerations let us now go to our first problem. The hamiltonian of eq. (||) can be decomposed 
as: 

H = Ho + V (18) 

V = W{r) (20) 

allowing the definition of the Lippmann-Schwinger basis. 

Its essential property is its spherical symmetry. We will see how this symmetry leads directly to the result above, 
avoiding the diagonalization procedure used in section II. In order to conserve this symmetry it is necessary that the 
nucleus prepares only spherically symmetric states. Thus, we can foresee that the CSCO corresponding to the pointer 
basis must contain the generator of angular rotation so it must be {H, L^, L^}. The pointer basis must then be the 
usual {jw,?,™ >} basis. 

In fact, let us consider an initial state functional po with spherical symmetry. If L is the generator the of densities 
rotations, a rotation of the state by an angle Ip gives po = exp{— i h -Ip} po, or equivalently 

(exp{-iL • ^}/9o |0) = {po\ e^p{i L ■ Tp} O exp{-i L ■ Ip}) , 

for any observables O G C In the last expression, L is the angular momentum operator. From the last equation we 
obtain 

(Lpo|0) = (po|[I,0])=0. (21) 

The observables O have the form 



= ^^ J dujOimd'm'{(^)\(^lm){ujl'm'\ + 

I m I' m' 



1^1^ I I dLjdLj'Oi^j,^'{uj,uj')\ujlm){uj'l'm'\, (22) 

/ m I' m' 



and therefore equation (^ gives 
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{po\[L,\ulm){uj'l'm'\]) =0. (23) 

These equations are equivalent to 

{po\[L,,\ujlm){oj'l'm'\])^0, {po\[L±,\ujlm){oj'l'm'\]) = 0, L±^L.,±Ly. (24) 
Taking into account that 



L±\l,m) = ^Jl{l + I) - mim ±l)\l,m±l) , (25) 

we obtain 

= (m — r7i')(po| 1^^, m) {uj' , I' , m'|) 



= + 1) - m(TO± 1) {po\\uj,l,m±l){uj\l' ,m'\) (26) 
-^l'{V + l)-m'{m'Tl) (polk, I, m){u;', I', m' T 1|). 
both for UJ = Lj' ov Lj ^ Lo' . These equations give 

{po\\uj,Um){uj' ,l\m'\) = p%uj,uo') Sw S^rw (27) 

(po||w,Z,TO)(w,/',m'|) = p^i{uj) 5iv Synm'- 

Thus any symmetric po is diagonal, in I and m. We will repeat the conclusion of the previous section and see, in this 
particular case, how the time evolution produces the diagonalization: namely p['(w,a;') will vanishes when t — s- oo. 
As in eq. (p^), the time evolution reads: 

= {pt\0) = {poW'^'Oe'^"') = (28) 

+1 



1=0 m=-l 

/r. OO +^ 

dbj / dcj'^ p°(w, w') exp{i(a; - tj')i} ^ 0;,„^;,„(a;, tj'). 



;=o m=-; 
Riemann-Lebesgue theorem can be used in this expression to obtain the "final" state 



(0)p. = (p.|0) = }^{pt\0) = j ckuf^p^iio) J2 Oir,.,iM. (29) 



1=0 m=-l 



where p^{uj,uj') has disappeared. If we define^ in analogy to eqs. (^ and (|7|) the functional {Lj,l,m,\ acting on an 
observable O of the form given in equation (E2h by (tj, I, m\0) = Oi„i^im{'^), we can give the following expression for 
the asymptotic form of the state 

/oo +/ 

1=0 m=-l 

Of course (p, | is also spherically symmetric since it is symmetric in / and m. We do not address in this paper the 
mechanism used by the nucleus to prepare the states in such a spherically symmetric way. But we are just studying the 
case where all the elements of the nucleus are spherically symmetric, and also the quantum states involved, because we 
are precisely trying to explain the breaking of this symmetry and the appearance of the radial structure. Then, it is 
clear that the only possibility is to begin with a spherically symmetric structure and with states satisfy the equations 
above. 

Therefore: 

i. -The origin of the decoherence in the energy is the time evolution. 

ii. - The origin of the decoherence in the angular variables is the preparation of the quantum state, which is spherically 
symmetric in our model. This symmetry is preserved by the time evolution. 

Let us observe that, as our model is spherically symmetric, any CSCO {-ff, L^,Lz}, for any arbitrary z axis, will 
correspond to the final pointer basis. But if the symmetry would be cylindrical along the axis z, being the center 
of an angular coordinate tp with generator L^, the only CSCO related with a final pointer basis, for the cylindrical 
quantum states, would be {H,pz, Lz}. So we see how the symmetry of the equation and the states defines the final 
pointer bases and their number. 
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B. The limit — > and the classical p)^''\q,p). 

Let us now compute the classical analogue of p, , as promised at the end of section II. We will prove that the 
distribution function pl'^' (g, p) , that corresponds to the density matrix via the Wigner integral [l^ is simply a 
function of the classical constant of the motion, in our case H{q,p), L^(q,p), Lz{q,p), precisely: 

p'^\q,p) = p.{H{q,p),L\q,p),L,{q,p)) (31) 
To simplify the demonstration let us consider only the constant H{q,p). From eq. ( p^ we have: 

p* = y pt:{uj){iLj\duj (32) 

So we must compute: 

pW(g,p)-7r-l f{u;\\q + X){q-X\)e^^P^dX (33) 



We know, from ||^, section II C, that the characteristic property of (w| is: 

(wli?") = u;" (34) 

Using the relation between quantum and classical inner products of operators ( eq. (2.13)) we deduce that the 
characteristic property of p^\q,p) is: 

p(f\q,p)[H{q,p)rdqdp - c." + 0{h) (35) 

for any natural number n. From now on we will work in the limit /i — > and therefore all the 0{h) will disappear. 
Thus Pu^\q,p) must be: 

p^f\q,p)=5{H{q,p)-u;)>0 (36) 

which turns out to be a distribution, namely a functional as p,. Therefore, going back to eq. (^^ and since the 
Wigner relation is linear, we have: 

P^^\q,p) = / p.iLu)p^J^Hq,p)dLU = / p,iLu)diHiq,p)-Lu)dLU - p.{H{q,p)) > (37) 



q.e.d. 

Generalizing this reasoning we can prove eq. (|3l|). Moreover the generalized eq. (|37| ) reads: 



iq,P) = J ^ P*{^J,'nfi)p^^i]^{q,p)duj (38) 

l,m 



where p,(w, I, m) = p'^{lu) of eq. (|2^) [] and pi^.l„(g,p) reads: 

pZUi^p) = S(H{q,p) - Lo)S{L'iq,p) - l{l + 1))5{LMp) - m) (39) 

and can be interpreted as the state where uj,l, m are well defined and the corresponding classical canonically conjugated 
variables completely undefined since p^^|„ is not a function of these variables. Then there is a complete coincidence 
with the result of the previous section. From the last two equations we obtain ( ^ as we promised. Now all the 



'^Even if, by symmetry reasons m is absent as an index of Pj'(w) we have introduced this index in p,(cj,/,m) for two reasons: 

i. -It may be present in a more general case (as the cosmological one). 

ii. - It is present in p^'^'/^. 
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classical canonically conjugated variables of the momenta H, L^, Lz do exist since they can be found solving the 
corresponding Poisson brackets differential equations 0. But as the momenta H,L^,Lz, that we will call generically 
I, are also constants of the motion, we have I' = —dH/da = 0, where a is the coordinate canonically conjugated to I, 
so H is just a function of I and: 

a* - = w{l) = const. (40) 

ol 

so: 

a = rzi{l)t + aa (41) 

These are the classical motions corresponding to the motion of the wave packet of the previous subsection. As in 
this section I is completely defined and ao completely undefined due to spherical symmetry, in such a way that the 
motions represented in the last equation homogeneously fill the surface I = const, (really H = const., = const., 
Lz = const, for our case), namely a usual classical torus of phase space. 

Then, eq. ( |3^ ) can be considered as the expansion of pi'^^\q,p) in the classical motion just described, contained in 

Pni mil-^P)^ each one with a probability p*(n, l,m). 
Summing up: 

i. - We have shown that the density matrix p[t) evolves to a diagonal density matrix p^, 

ii. - This density matrix has pi'^^\q,p) as its corresponding classical density. 

iii. - This classical density can be decomposed in classical motions where H, L^, and Lz remains constant, answering 
question i of the introduction. 

iv. - Finally, for r ^ oo the classical analogue of eq. (p^): 

show that in the classical motions we can consider = when r ^ oo |^. Thus the classical motions far from the 
nucleus can be considered as radial, answering question ii of the introduction. 
Finally as: 

[9 - 6i("))(5((p - ip'-°^)d9'-°Uip^°^ = 1 (43) 
where 9 and ip are the usual polar coordinates. Since far from the nucleus Z = m = 0, we can write eq. 

pl^'(g,p)= f p,{iu,0,0)6{H{q,p) - iu)6{9 - 9(°^)6{^ ~ ^^°^)d9^°^d^(''Uio (44) 



where 5{H{q,p) - uj)6{9 ~ 9'-°^)6{ip - ^(o)) can be considered as the classical density of the classical particles with 
radial motions defined by the energy uj and the angles 6''^*'^ and (p^^\ The final classical equilibrium density is thus 
decomposed in the densities of radial classical trajectories, as in eq. (||). 

So finally we have an isotropic ensemble of particles in radial motion. This is the classical statistical limit of our 
wave function. We have gone from quantum mechanics to classical statistical mechanics [^6|. To single out any one 
of the trajectories to obtain a single classical motion is clearly impossible, at the statistical mechanics level, since it 
would brake the spherical symmetry of the initial wave function. But the statistical state is composed of single radial 
motions, at the classical level, as any population is composed of individuals. 



C. Localization and correlations. 

These phenomena will appears for adequate potential and initial conditions (see [Q). U W{r) = we know that 
the wave packet will spread, so the real W{r) must be such that the eventual spreading be as slow as necessary in 
order to see the radial trajectories in the bubble chamber. 



*On the contrary, H and have not quantum canonically conjugated dynamical variables since their spectra are bounded 
for below. 

®A close study of the initial conditions, that necessarily are located near the center of the nucleus will improve the demon- 
stration of this result. For the sake of conciseness we do not include these reasonings. 
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D. Discussions and comments. 



1. The time. 



We have postulate the global existence of time in section IIA. This postulate was motivated in the fact that this is 
a global feature of the universe. But, in order to mimic the cosmological case of the second example, let us imagine 
for a moment that our model would be considered as a model of the whole universe Then it would be impossible 
to postulate the existence of time based in an exterior global structure. A way to solve this problem would be to 
postulate decoherence through eq. (p^, namely the existence of some kind of evolution with such a final diagonal 
state which corresponds to the fact that the universe really ends in a classical states. Then we can define the time 
over the classical trajectory labelled by uj and I via eq. (E2|) as: 



M- 



'di 



I2M 



2Mr^ 



Wir) 



(45) 



1. e.: 



1 



dr 



1(1+1) 



W{r) 



(46) 



Thus we can see the interrelation of the time, the decoherence, and the elimination of the uncertainty relations. 
The usual procedure is to postulate the existence of time and follow the chain time— >decoherence^elimination of the 
uncertainty relations. In this case we are using the traditional way of thinking of classical and non relativistic quantum 
mechanics: time is a primitive concept. But, as we have explained, another chain is feasible: decoherence— ^elimination 
of the uncertainty relations— >time, and that would be the chain that we must use if we consider the whole universe. In 
this last case we could postulate (based on a obvious observational fact) that the quantum physics of the universe is 
such that it has a tendency towards a classical regime (decoherence + elimination of the uncertainty relations). This 
tendency would be the primitive concept in this case and time would be a derived concept in perfect accord with Mach 
philosophy (see also |^l|). So question iii of the introduction is answered. We will come back to these arguments 
in the cosmological case in section IVD3 where we will go further on. We will postulate the existence of a parameter 



ry, that at the quantum level will take the role of time, namely \ri) 
this parameter that, therefore, becomes a classical one. 



iHt 



|0), and show that there is decoherence in 



2. Why the outgoing solutions? 

The answer to this question can be searched in one of the essential properties of the actual state of the real universe: 
it is time- asymmetric, in spite of the fact that its evolution laws are time-symmetric, e. g.: 

i. - Even if we have incoming and outgoing scattering states, only the second ones behave spontaneously, while the 
first ones must be produced by a source of energy 

ii. - Even if we have advanced and retarded solutions of the Maxwell equations the time asymmetry of the state of 
the universe force us to use the later ones and to neglect the former ones. 

iii. -At the quantum level, the space of admissible solutions is not the whole Hilbert Space but a subspace of this 
space, where causality can be introduced Fn. 



This will be the case in the first problem. In the second one we will deal with the real universe and reproduce the same 
kind of arguments. 

'^^Incoming states will transform the stable nucleus in a unstable one. Outgoing states will radiate the energy produced when 
the unstable nucleus evolves to a equilibrium state in an spontaneous way. 

^^As causality is related with the analytic properties of wave functions , when we promote the energy lj to a complex variable 
z, we can choose the admissible solutions as those wave functions that are analytic and bounded in the lower complex half-plane 
of variable z, even if we could choose those that have the same properties in the upper plane(or more precisely such that they 
belong to the Hardy class function space from below (or above) PTl)- ^^i this way causality is introduced and its physical 
consequences: dispersion relation, the fluctuation-dissipation theorem, the growing of entropy, etc. 13. 
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Namely, as the laws of physics are time symmetric they always give us two t-symmetric possibilities. The universe 
is t-asymmetric and corresponds to one of this possibilities. 

From what we have said it must be clear that, in order to be completely satisfactory, the choice must be done in 
the whole universe, namely it must be global [2^. Let us only consider the case "i", we must only choose outgoing 
states answering question iv. 

For our first problem we may use the WKB solution of the equation: 



+ k^{x)i){x) =0 



which is: 

7/1(2;) = [fc(a;)]~^ exp[±i j k{x)dx\ 
Then if we write the solutions of the eq. (0) as: 



the function u^i{r) satisfies the equation: 

d'^u^iir) 



dr^ 



2M 



uj ' ^ „ - W(r) 



Uujiir) = 



If we can call: 



the WKB solution is 



kliir) = 2M 



2Mr^ 



W{r) 



Ar) ^ [kuji{r)] 2 exp[±i / k^i{r)dr 



(47) 



(48) 



(49) 



(50) 



(51) 



(52) 



These are the two possibilities produced by our time symmetric theory. Then, in order that solution ( p2[) would 
be out-going we must choose the sign +. In fact, far from the nucleus, when t —>■ oo,r ^ oo we have kuii{r) — > fc^ = 
+\/2Muj and considering the time evolution factor we have: 

e~''"^u^i{r) = [kujY^ exp[-i(a;t - k^r + cmst.)] (53) 
and the wave evolution becomes the outgoing unilateral shift: 



r = - — t + const. 

Kin 



(54) 



So the question iv of the introduction is also answered |^. 



Perhaps it is interesting to observe that, following the ideas of paper ^^li this typical outgoing nature of the shift (that is 
extremely important in the cosmological models |^^) it is only possible because, from the very beginning the configuration space 
has a characteristic structure: it is spherically symmetric and this fact defines the asymmetry r > 0, namely the asymmetry that 
says that the origin O is substantially different than the sphere at the infinity. Therefore; it is an asymmetry in configuration 
space that really introduces time- asymmetry j2^. Without this asymmetry all the treatment would be impossible, since the 
"outgoing" notion itself would be meaningless, and the choice of the lower half-plane unmotivated, (at least in the case where 
we consider our system as the whole universe and therefore we have no other physical phenomena to play with). Therefore the 
fact that an isolated nucleus radiates and never receives spontaneously radiation from the exterior, in a conspirative way, is 
a consequence of the fundamental time-asymmetry of the universe. But if, as a theoretical example, we consider this isolated 
system as the whole universe it is a consequence of the asymmetry of the configuration space of the model, in complete 
agreement with ref. [Q. In the second example the time asymmetry has, more or less, the same origin; the global asymmetry 
of the universe as explained in reference [M. 
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IV. THE COSMOLOGICAL PROBLEM. 



A. The model. 

Let us consider the flat Roberson- Walker universe ( ||2^, [Q) with a metric: 

ds^ =a^{r]){dr/^ -dx^ -dy^ -dz^) (55) 

where rj is the conformal time and a the scale of the universe. Let us consider a free neutral scalar field $ and let us 
couple this field with the metric, with a conformal coupling = i). The total action reads S =^ Sg + Sf +Si and the 
gravitational action is: 

Sg^M^ J dr^[-^a -V{a)] (56) 

where M is the Planck mass, a= da/drj, and the potential V contains the cosmological constant term and eventually 
the contribution of some form of classical mater. We suppose that V has a bounded support < a < ai. We expand 
the field $ as: 

a>(,y,x) = j he-'^'-^dk (57) 

where the components of k are three continuous variables. 

The Wheeler De-Witt equation for this model reads (compare with eq. (|l|) for the first example)): 

H^{a,<i>) = {hg + hf + hi)^{a,<i>) = (58) 

where: 

hs = ^dl + M'V{a) 
hf = -\ j{dl-efl)dk 

1 ,2„2 / j-2. 



h, = -m'a' j f^dk (59) 

being m the mass of the scalar field, k/a is the linear momentum of the field, and 5k =9/i9/k- 
We can now go to the semiclassical regime using the WKB method ( (25)), writing ^'(a, $) as: 

^{a, $) = exp[iM^S{a)]x{a, $) (60) 

and expanding S and x as: 

S = So + M-^Si + x^Xo + M-\i + ... (61) 



To satisfy eq. (58) at the order M , S{a), the principal Jacobi function, must satisfy the Hamilton- Jacobi equation: 

2 



dS\ 
da J 



= 2V{a) (62) 



We can now define the (semi) classical time as in section III C (but now using an approximate solution only). It is a 
parameter rj = ri{a) such that: 

d dS d 



drj da da ^ ^^da ^ ^ 



So: 
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da (64) 



V2Jao VW) 

that must be compared with eq. (^), but m this equation the trajectories are completely classic, while in eq. ( |64 
only a is classic while $ remains as a quantum variable. The solution of eq. ( |64|) is a = iF{ri, C), where C is an 
arbitrary integration constant. Different values of this constant and of the ± sign give different classical solutions for 
the geometry. 

Then, in the next order of the WKB expansion, the Schroedinger equation reads: 

= Hv)x (65) 

where: 

h{r)) ^ hf + hi{a) (66) 

precisely: 

92 



dk (67) 



where: 

= rn^a^ + fc^ = rn^a^ + w (68) 

and vj = k"^ and k — |k|. So the time dependence of the hamiltonian comes from the function a = &(??). 

Let us now consider a scale of the universe such that aout 3> ai. We will consider the evolution in this region 
where the geometry is almost constant. Therefore we have an adiabatic final vacuum |0) and adiabatic creation and 
annihilation operators and at- Then h = h{aout) reads: 

^^al^ai^dk (69) 

We can now consider the Fock space and a basis of vectors: 

|ki,k2,...,k„,...)-|{/c})=a|,^a[^...a[^^...|0) (70) 
where we have called {k} to the set ki, k2, k„, ... The vectors of this basis are eigenvectors of h : 

h\{k}) ^u;\{k}) (71) 

where: 

io= ^-^^ E (72) 
ke{k} ke{k} 

We can now use this energy to label the eigenvector as: 

|{fc}>-k[k]) (73) 

where [k] is the remaining set of labels necessary to define the vector unambiguously, {jcj, [k])} is obviously an 
orthonormal basis so eq. (69) reads: 



io\uj,[k]){uj,[k\\dujd[k] (74) 
This is the hamiltonian that corresponds to (^ in the cosmological case. 
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B. Decoherence in the other dynamical variables. 



We can obtain the decoherence in the energy as in section II. Then, if we reintroduce the other dynamical variables 
in eq. (113) we obtain: 



{pA = j Pc.[k][k'](^, [k], [k']|dc^d[k]d[k'] (75) 

where {((^j fk , [k']|, {uj,uj', [k], [k']} is the cobasis {(wj, (w, w'j} but now showing the hidden [k]. This equation corre- 
sponds to (|l4) in the cosmological case. 

Let us observe that if we would use polar coordinates for k eq.(|5^) reads: 

<^{x,n) = J ^(jiumdk (76) 

Im 

where: 

<j^ki^ = fkAv,r)Y^{e,^) (77) 

where fc is a continuous variable, / = 0, 1, ; m = — I] and Y are spherical harmonic functions. So the indices 
k, 1,171 contained in the symbol k are partially discrete and partially continuous. 

As pi — p, then p*[k'][ki — Puj[k;][k;'] s-nd therefore a set of vectors {|w, [1])} exists such that: 

Pc^[k][k']k, [l])[k']rf[k'] = [l])[k] (78) 

namely [1])} is the eigenbasis of the operator P(j[k][k']- Then p^^jij can be considered as an ordinary diagonal matrix 
in the discrete indices like the / and the m, and a generalized diagonal matrix in the continuous indices like k 
Under the diagonalization process eq. ( [TS] ) is written as: 

(p*i = / c^[kiv.[i][i']<^]f/[i!!] '(t^, [1"], [i'ii<;'idu;d[k]d[k']d[i]d[iv[i'']dn (79) 

where f7j^|' is the unitary matrix used to perform the diagonalization and: 



where: 



so we can define: 



= Puj[\\5[\][i'] (80) 
= P^[\] = / t^mV^[k][k']f^m'''rf[k]rf[k'] (81) 



(a;, [1]| ^ {u;, [1], = J [/^(c., [k], [k'] |[/ffl^d[k]d[k'] (82) 

We can repeat the procedure with vectors (w, w', [k], [k']| and obtain vector (tjjtj', In this way wc obtain a 
diagonalized cobasis {(cj, (uj,uj', [1]}. So we can now write the equilibrium state as: 

P* = / P^u]{uJ,m^d[l] (83) 



^''E. g.: We can deal with this generalized matrix rigging the space S and using the Gel'fand-Maurin theorem this 
procedure allows us to define a generalized state eigenbasis for system with continuous spectrum. It has been used to diagonalize 



hamiltonians with continuous spectra in Q, etc. 
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which corresponds to ( p^ ) in the cosmological case. Since vectors (w, [1]| can be considered as diagonals in all the 
variables we have obtained decoherence in all the dynamical variables. This fact will become clearer once we study 
the observables related with this vector and introduce the notion of final pointer basis. 

So, let us now consider the observable basis {\u!, [1]), [1])} dual to the state cobasis {{uj, {uj,uj', From 

eq. and as the u> does not play any role in the diagonalization procedure we obtain: 

= |a;,a-',[l])=K[l])(c.',[l]| (84) 

So in the basis {\u!, [1]), \u!,u!', [1])} the hamiltonian reads: 

h = J uj\uj,[l])dLud[l]^ J uj\uj,[l]){uj,[l]\dujd[l] (85) 

Now, we can also define the operators: 

L = J l\uj,[l])dujd[l]= J l\uj,[l]){tJ,[l]\dujd[l] (86) 



that can also be written, as in eq. ([171): 

L, = fk\co, [l])dLud[l] = / k\LU, [1]){lu, [\]\dujd[l] (87) 



where i is an index such that it covers all the dimension of the 1 f^. Now we can consider the set {h,Li), which is 
a CSCO, since all the members of the set commute, because they share a common basis and find the corresponding 
eigenbasis of the set, precise \uj, [1]) since 0: 

h\uj,[\])^i,\uj,[l]) 



L,\cj,[l])^l,\u;,[l]) (89) 

Of course the Li are constant of the motion because they commute with h. From all these equations we can say that: 

i. - {h, Li) is the final pointer CSCO. 

ii. - {\ui, [1]), \uj,uj', [1])} is the final pointer observable basis. 

iii. -{(w, (w, w', [1]J| is the final pointer states cobasis. 

In fact, from eq. (p3) we see that the final equilibrium state has only diagonal terms in this state (those corre- 
sponding to vectors (w, [1]|) , it has not off-diagonal terms (those corresponding to vectors {lu,uj', [1]|,(w, [k], [k']|, or 
{u!,u!', [k], [k']|), and therefore we have decoherence in all the dynamical variables. 



C. The limit ^ and the classical pi'^'\x,k). 

Let us restore the notation {1} — {u;, [1]), {k} — (uj, [k]), as in eq. ( [73| ) and let us consider the configuration kets 
\{x}) — \r], [x]). Since we are considering the period when a ^ Uout the system with hamiltonian (|6^) is just a set of 
infinite oscillators with constants (flout) that represent a scalar field with mass maout- Then we are just dealing 
with a classical set of N particles, with coordinates [x] and momenta [k] . Then, as in eq. (^) we can introduce the 
Wigner function corresponding to generalized state |{/}). 

p;^)([x], [k]) = ^-4^ j ({/}|x + A)(x - A|)e2*W-Md4"A (90) 

Using the same reasoning that we have used to obtain eq. (pfl) it reads: 



In principle the matter field "I> may have any number of particles A*'. But since we are working in the final stage of the 
universe evolution with a ~ aout, this number is a constant. Then the number of observables in the CSCO is 4A'' and the ket 
in configuration variables read \rj, [x]) = [{i}), where [x] is the space position of the AN particles. 
^^In some occasions we will call h = Lq and lo = Iq. 



14 



p;^)([x],[k]) = n'^(ir([x],[k])-y (91) 

i 

where Lj^([x],[k]) is the classical observable obtained from Li via the Wigner integral (considering h = Lq and 
including among the indices i). Now, with the new notation (pO|), eq. (pi) reads: 



P*{l}m\d{l} (92) 



pl'^^([x],[k]) =7r-4^y'({p,|x + A)(x- A|)e2'W-Md4"A (93) 



^^^([x], [k]) = pr^(C([x], [k]),Lr([x], [k]), ...) (94) 

So finally: 

^(^)([x],[k])~ I d{Z}p{,jpr)([x],[k])5({L^}-{Z}) = I d{?}p{,|in'5(Lr-^.) (95) 



the if we call: 



we obtain, as in eq. (Bl^ 



P* 



The last equation can be interpreted as follows: 

i. - S{{p} — {I}) is a classical density function, strongly peaked at certain values of the constants of motion {I}, 
corresponding to a set of trajectories, where the momenta are equal to the eigenvalues of eqs. (|8|) and (p9|), namely 
L^ = h (z = 0,1,2,...). 

ii. - is the probability to be in one of these sets of trajectories labelled by {I}. Precisely: if some initial density 
matrix is given, from eq. ( |83| ) it is evident that its diagonal terms p^iy are the probabilities to be in the states (w, [1] | 

and therefore the probability to find, in the corresponding classical equilibrium density function pl^''([x], [k]), the 
density function SijL^'} — {I}), namely the probability of the set of trajectories labelled by {/} = (w, [1]). 

iii. - As in eq. ([40[ ) let a be the coordinate classically conjugated to 1 and let be Hq the coordinate a at time 7] = 0, 
then we obtain the classical trajectories: 

a = l?7+ao (96) 

iv. - Let us now call p*{l} = P{i}[a.o] - Really P{/}[ao] is not a function of ap, it simply is a constant in ag, since ao is 
only an arbitrary point and our model is spatially homogenous. Then we can write: 



P{i}[a 



in this way we have changed the role of slq , it was a fixed (but arbitrary) point and it is now a variable that moves all 
over the space. Then eq. (p5|) reads: 



So if we call 



we have: 



pi''\[^], [k]) ~ / Piin.o]I[^iL7 - h) n S{aj - aoM^oW} (98) 



p;fj)^^,([x], [k]) = \{6{LY k) n 5{a, - ao,) (99) 

1=0 j = l 



pl^')([x], [k]) ~ / P{,}[aolpSU](W. [k])rf[ao]d{0 (100) 



From eq. ( p9[ ) we see that P|'Jjj^^j([x], [k]) 7^ only in a narrow strip around the classical trajectory ( p6[ ) defined by 
the momenta {/} and passing through the point [ao] (really the density function is as peaked as it is allowed by the 
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uncertainty principle, so its width is essentially a 0(h ), since the 5— functions of all the equation are really Dirac's 
deltas only when h ^ 0) . So we have proved eq. ( |lOO| ) which, in fact, it is eq. as announced 

Then we have obtained the classical limit. When ?7 — > oo the quantum density p becomes a diagonal density matrix 
/9» . The corresponding classical distribution pi'^^^ ( [x] , [k] ) can be expanded as a sum of classical trajectories density 
functions /9^'j'j''j^^j([x], [k]), each one weighted by its corresponding probability P{;}[ao]- So, as the limit of our quantum 
model we have obtained a statistical classical mechanical model, and the classical realm appears. 



D. Localization and correlations. 



Under adequate initial conditions the motion can be concentrated in just one trajectory showing, the presence 
of correlations in this trajectory Q. The evolution of the concentration depends on potential V{a). Of course, our 
"trajectories" are not only one trajectory for a one particle state, but they are N trajectories (each one corresponding 
to a momenta (Zi, I2, ■■■In) = {1} and passing by a point (ai, 3.2, a„) = [a]) for the n particle states. 



E. Discussion and comments. 



!■ Characteristic times ■ 



The decaying term of eq. pO) (i- e. the second term of the r. h. s.) can be analytically continued using the 
techniques explained in papers [||, jl^, and [Q. In these papers it is shown that each pole Zi = uii — iji, of the 
S-matrix (corresponding to the evolution ai„ aout, see |2^) of the problem considered, originates a damping factor 
e Then if 7 = min(7i) the characteristic decoherence time is 7^^. This computation is done in the specific models 



of papers |2j]. 
since 7"-^ » 1. 



If 7 <C 1, even if the Riemann-Lebesgue theorem is always valid, there is no practical decoherence 



2. Sets of trajectories decoherence. 



It is usual to say that in the classical regime there is decoherence of the set trajectories labelled by the constant of 
the motion uj, [1]. This result can easily be obtained with our method in the following way. 

i.- Let us consider two different states and jij^'il']) that will define classes of trajectories with different constants 
of the motion {u, [1]) 7^ (w', [!']). We must compute: 



{Lj[i]\p,\ij'[i']) = ip^Wivcj'm']) = 



]{uj"[l"]\duj"d[l" 



/[!][!']) = 



(101) 



due to the orthogonality of the basis {(oj, [1]|, (aj,aj', [1]|} . 
ii.- But if we compute: 



p^„[y,]{uj"[l"]\duj"d[l" 



41]) = 



^''in this section, as in sectionn IIB, we have faced the following problem: 

pl'^''([x], [k]) is a a constant that we want to decompose in functions pj'j'jj^^j ([x], [k]) which are different from zero only around 
the trajectory (^6|) and therefore are variables in a. 
Then, essentially we use the fact that if f{x,y) — g{y) is a constant function in x we can decompose it as: 



9iy) = J 9{y)5{x - xo)dxo 

namely the densities 5{x — xq) are peaked in the trajectories x = xo = const^,y = var. and, therefore, are functions of a;. This 
trajectories play the role of those of eq. (p7|). 

As all the physics, including the correlations, is already contained in eq. (p5|), the reader may just consider the final part of 
this section, from eq. (OT) to eq. ( IOC ) a didactical trick. 
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J P."li"]Siu; - cj")Si[l] - [l"])cL;"d[l"] = ^ (102) 

The last two equations complete the demonstration. We will discuss the problem of the decoherence of two trajectories, 
with the same {1} but different [ao] in subsection 1VD4. 



3. A discussion on time decoherence. 



It is well known that one of the main problems of quantum gravity is the problem of the time definition (see p9| ) . 
A not well studied feature of this problem is that, there must be a decoherence process related with time, since time 
is as a classical variable. In this subsection, using the functional technique, we will give a model that shows that this 
is the case (but we must emphasize that this subject is not completely developed). 

Let us postulate that there is a parameter rj such that the quantum states evolve as 0: 

\V) = e-'^'^m (103) 

We must compute (?7|p*|?7') where |r/) and I77') are two states of the system for different times. We do not know if 
{''l\p*W) "^ill decohere or not. If it decohers we can say that parameter 77 is a classical one. \'ri){ri'\ can be considered 
as an observable, then: 

{V'\P*\V)^{P*\\V)W\) (104) 

But: 

{LuMirj'D (c.|e-''''|0)(0|e""'') = [e''"''(c.|e-^'']||0)(0|) (105) 
Now, for any observable O we have: 



[g,M'(^|e-«M]|[yo^,|^')rf^' + ... ^ {uj\[J 0^,e~"^''^uj')e"^''i'duj']) = 

where the second term disappears since {uj\lu',uj") — 0. Thus: 

(c.|h)(7y'|) = e— ("'-")(c.||0)(0|) (107) 

So now we can compute the following two cases: 
i.- 

{V'\P*\V) = iP*\\v){v'\) = [ f pU^\dLo]\\r,){r,'\) = f p^e"- (c.| |0) (0|)d^ ^ (108) 



when 1 77' — r/l — !■ 00 , due to the Riemann-Lebesgue theorem, 
ii.- Analogously: 



{V\P*\V) = J P.(^||0)(0|)dc.^O (109) 
So we have time decoherence for two times 77 and rj' if they are far enough. 
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Of course 77 is the conformal time of eq. (pq), since (103) is a consequence of (p5|). But now we have postulated this last 



equation and we are searching the quantum properties of 77. 
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This result is important for the problem of time definition, since in order to have a reasonable classical time this 
variable must first decohere. The result above shows that this is the case for 77 and rj far enough 0, but also that, for 
closer times (namely such that their difference is smaller than Planck's time) there is no decoherence and time cannot 
be considered as a classical variable. Classical time is a familiar concept but the real nature of the non-decohered 
quantum time is opened to discussion. But we must remark that, some how, we have followed the second line of 
thought of section llCl: wc have suppose the existence of an evolution e"*'"', where r] is only a parameter. We have 
prove that decoherence appears and find that 77 behaves like a classical variable. May be this is the better way to 
introduce the classical time to postulate a "quantum" t] and to find its properties. Moreover we have proved that the 
second line of section IICl can as well be followed. 

J^. Decoherence in the space variables. 

Now that wc know tha t the re is time decoherence we can repeat the reasoning for the rest of the variables a at time 
77 = and changing cq. ( |l03[) by: 

|[a]) =e'W'[il|0) (110) 

and we will reach to the conclusions: 
i.- 



([a]|p*|[a'])^0 (111) 



when |a — a'| — > c». 
ii,- 



([a]|p*|[a])^0 (112) 

therefore there is also decoherence between two trajectories with the same {1} but different [ao]. 

These facts complete the scenario about decoherence and the final classical limit. An analogy of section IIIC2 in 
the cosmological case can be found in paper |pT[| . 

V. CONCLUSION. 

We are convinced that the method of papers [|| and |l^ is the best way to study both the Mott and the cosmological 
problem and to find the analogies and differences between then. We hope that the reader will share our conviction. 

For the most important of both model, the cosmological one we have shown that after the WKB expansion and the 
decoherence and the final classical limit process our quantum model has: 

i. -A defined classical time 77 and a defined classical geometry related by eq. (|65|). 

ii. - Decoherence has appeared in a well defined final pointer basis. 



iii.- The quantum field has originated a classical final distribution function (eq. (100)) that is a weighted average 
of some set densities, each one related to a classical trajectory. The weight coefficients are the probabilities of each 
trajectory. 

We can foresee that if instead of a spinless field we would coupled the geometry with a spin 2 metric fiuctuation 
field the result would be more or less the same. Then the corresponding quantum fluctuations would become classical 
fluctuations that would correspond to matter inhomogeneities (galaxies, clusters of galaxies, etc.) that will move 
along the trajectories described above. But this subject will be treated elsewhere with greater detail. 



'Using the method of section IVCl we can compute 7. Decoherence will take place for \ri ~ r]'\ > 7 



-1 
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